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We present a combined QCD analysis of diffractive and leading proton deep inelastic scattering data using the framework of fracture 
functions. It is shown that this framework allows a precise and unified perturbative QCD description for the data, alternative to 
those that relay on Regge factorization. 


1 Introduction 

In recent years considerable attention has been paid to 
diffractiv^ jEiid leading-proton deep inelastic scattering 
processesld’Bfllj. These two kinds of processes differ con¬ 
siderably in the kinematical regions where they are pro¬ 
duced and also in the non-perturbative mechanisms or 
models used to explain them and, consequently, their 
description is usually done within completely unrelated 
frameworks and in terms of different structure functions: 
the leading-proton and the diffractive deep inelastic scat¬ 
tering structure functions, and respec¬ 

tively. 

However, if the experiment under consideration is 
able to identify the final-state proton in the diffractive 
processes, as in the ZEUS measurements, or if there is 
confidence in the dominance of the single dissociative 
process 7 *p ^ Xp, both kinds of events can be thought 
as of semi-inclusive nature, with identical final-state par¬ 
ticles produced in the target fragmentation region. From 
this point of view, the perturbative QCD framework for 
their description must be identical, with only one factor¬ 
ized observable, even though the specific models for their 
non-perturbative features are completely different. 

In perturbative QCD, the most appropriate descrip¬ 
tion, for semi-inclusive DIS events in which the identified 
final-state hadron is produced in the target fragmenta- 

t on region, is the one that includes fracture functions 
These functions can be understood essentially as par- 
ton densities in an already fragmented target, and extend 
the more familiar description of semi-inclusive events in 
terms of parton densities and fragmentation functions, 
allowing a leading-order description for target fragmen¬ 
tation events in the forward direction, and solving some 
problems asociated to factorizationIjCl. At the same time, 
the formalism provides not only the motivation for es¬ 
tablishing a connection between the kinematical regions 


of these kinds of processes, but also rigorous predic¬ 
tions about the perturbative QCD behaviour of the cross- 
sections. A priori this behaviour is not the same as that 
of the ordinary structure functions due to appearence of 
an inhomogeneous ter^in the Altarelli-Parisi equations 
for fracture functions ffl. 

In the following we present results from a QCD global 
analysis of recent semi-inclusive deep inelastic scatter¬ 
ing data produced by HI using the framework of frac¬ 
ture functions. It is shown that this approach unifies the 
description of diffractive and leading-proton phenomena 
and allows a perturbative QCD description without the 
usual assumptions about approximate Regge factoriza¬ 
tion. The resulting parametrization is also used to com¬ 
pute other observables measured by ZEUS, not included 
in the fit to HI data, finding also an outstanding agree¬ 
ment with the data. 


2 Definitions 


It is customary to define the leading-proton structure 

r p/Q\ 

function F 2 ^ ' from the corresponding triple-differential 
deep inelastic scattering cross section: 
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and the usual kinematical variables. 

Even though the processes accounted for in Eq. (1) 
are explicitly of a semi-inclusive nature, the formulation 
based on the leading-proton structure function is used 
instead of the usual approach for semi-inclusive deep in¬ 
elastic scattering in terms of parton distributions and 
fragmentation functions, because the last one only takes 
into account hadrons produced in the current fragmen¬ 
tation region and thus not contributing to the forward 
leading hadron observables. 

These problems, and those related to factorization 
of collinear singularities at higher orders, are overcome, 
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however, if the complete perturbative framework for 
semi-inclusive processes is taken into account, for which 
the LO expression for the production of very forward 
hadrons is 


j3 p 
“ ^target 

dx dQ'^ dz 


xQ'^ 



ejxMf^^{x, z, Q^) ( 2 ) 


where Mf^^{x,z,Q‘^) is the fracture function that ac¬ 
counts for target fragmentation processes and we define 
the variable z = Ep/ jEp as the ratio between the ener¬ 
gies of the final-state proton and the proton beam in the 
centre of mass of the virtual photon-proton system. For 
very forward protons then, 1 — ^ ~ z. 

Fracture functions can be thought of in terms of 
the elements of any lowest-order picture for hadroniza- 
tion, for example as the product of a flux of exchanged 
‘reggeons’ times their structure functions, or more for¬ 
mally as an ingredient of the perturbative QCD treat¬ 
ment: the non-perturbative parton distributions of a pro¬ 
ton fragmented into a proton. The latter choice can 
be generalized to higher orders, allowing a consistent 
analysis of the scale dependence of the semi-inclusive 
cross section. This scale dependence is driven by inho¬ 
mogeneous Altarelli-Parisi equations, reflecting the fact 
that the evolution is not only driven by the emission of 
collinear partons from those found in the target (homo¬ 
geneous evolution), but also by the fragmentation of par- 
tons radiated from the one struck by the virtual probe 
(inhomogeneous evolution). 

Defining the equivalent to Ei for fracture functions, 

= xY^e^,Mf^^{x,z,Q‘^), (3) 

i 

and taking into account the shift from z to the rela¬ 
tion between this function and the leading-raoton struc¬ 
ture function is quite apparent. In referenced, the use of 
fracture functions for the description of leading-hadrons 
produced in the target fragmentation region has already 
been discussed in relation to the analysis of very forward 
neutrons observed by the ZEUS collaboration. There 
it was shown that, neglecting contributions beyond LO 
coming from the current fragmentation region, what is 
usually defined as the leading-neutron structure function 
is just the fracture function contribution to the 
semi-inclusive cross section. 

Analogously to Eq. (1), the differential cross section 
for diffractive deep inelastic scattering is usually written 
in terms of the diffractive structure function E^^^^ 




dfd dQ^ dx^ 
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dissociative process 7 *p ^ Xp, implying that there is a 
proton in the final state, the contributions to it are again 
given by the same fracture function in (2), even though 
in a completely different kinematical region. 

Different kinematical regions correspond to different 
behaviours and also to different underlying models. The 
leading-proton struc^re function has been measured by 
the HI collaborationa, and has also been compared with 
predictions of different mechanisms, such as meson ex¬ 
change and soft colour interaction models, implemented 
in event generator programmes, none of which-, repro¬ 
duced integrally the main features of the data cl. The 
standard interpretation for the diffractive cross section 
is given in terms of ‘pomeron’ exchange. In this frame¬ 
work, different model estimates, and even QCD-inspired 
parametrizatipns of the ‘pomeron’ content, have been 
proposed however the comparison between model 
predictions and data again has been found to be rather 
poor unless a large number of additional elements is in- 
cludedU. 

In the language of fracture functions, both the 
leading-proton and the diffractive regimes are comple¬ 
mentary features of a more general semi-inclusive pro¬ 
cess. The approach, then, suggests and provides a bridge 
between the two regimes, which is particularly appropri¬ 
ate and even necessary, at least in the kinematical re¬ 
gion where neither the ‘pomeron’ nor the ‘reggeon’ ex¬ 
change picture alone are expected to describe the data, 
say Xjp ^ 0.1. 

3 Parametrization 

In order to obtain a parametrization for the proton-to- 
proton fracture function {(3,Q'^,Xjp) at a given ini¬ 
tial scale Qq, we select in the first place, a relatively 
simple functional dependence in the variables (3 and Xjp. 
If one were only interested in the diffractive regime, the 
natural choice would be a simple ‘pomeronic’ flux in Xjp 
times an ordinary parton distribution in /3. For the 
leading-proton regime the natural choice would be al¬ 
most the same, but with a standard meson or ‘reggeon’ 
flux or even better, something combining their effects. 
These kinds of parametrizations give relatively good ini¬ 
tial approximations to the description of the correspond¬ 
ing data sets; however their survival seems unlikely in a 
more precise analysis. 

In order to take into account small departures from 
the initial approximations, and also combine the two 
behaviours in such a way that for low Xjp (diffractive 
regime) the ‘pomeron’ picture dominates, while for low 
/3 and large Xjp the meson or ‘reggeon’ exchange picture 
emerges, we propose a modified flux such that the light- 
quark singlet component {Mq^^ = 3Mu^^ = = 


where Xjp = and the variable (3 is used instead of x. 
When this last cross section is dominated by the single 
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Figure 1: HI diffractive data against the outcome of the fracture 
function parametrization. 


Figure 2: HI leading-proton data against the outcome of the 
fracture function parametrization. 


Table 1; Parameters for Qq = 2.5 GeV^. 


Set 

01 jp 

Clp 

Plp 

ILP 

app 

Ng 

N, 

ag 

A 

-1.260 

14.395 

32.901 

2.627 

12.320 

0.041 

0.354 

0.450 

B 

-1.257 

12.556 

32.412 

2.338 

11.412 

0.047 

0.694 

0.648 


3Ms^^) of the fracture function is parametrized as 
xMP/P{P,Qlx^) = X (5) 

{C^ Pxy + Clp (1 - pypp (1 + aLp{l - x^Y^n] 

and similarly for gluons with the corresponding parame¬ 
ters Ng, Ug and bg. 

Even though at the initial scale Qq = 2.5 GeV^ the 
parametrization implies some sort of flux factorization, 
beyond the initial scale, the evolution equations drive 
the fracture function as a whole making the usual dis¬ 
crimination between ‘fluxes’ and ‘parton densities of the 
exchanged object’ somewhat ambiguous. 

In order to ensure that the evolution code is working 
properly in the region of large /3, we imposed a lower 
constraint on the (1 — /3) exponents for both quark and 
gluon distributions, 6^, bg > 0.1. Furthermore, in the 
expectation of very hard distributions we saturate this 
constraint and fix = 0, leaving only 8 free parameters. 
Doing this we obtain a global fit with XjQj^j/d.o.f.=1.09 
(y^ = 292.23, data points = 274), which we designate as 
Set A in Table 0. 

Similar parametrizations, but with softer gluons, 
yield slightly higher values, for example in Set B, 


where bg is set to 0.7 and bg remains 0.1, finding 
xLaz/d-0.f.=1.13. 

As is shown by the solid lines in Figures. 1 and 2, 
the accuracy of the fit is remarkably good in the case of 
HI diffractive data (x|fi£)/dato = 215.63/226), and also 
good for the leading-proton data. The dashed line in fig¬ 
ure 1 shows the contribution coming from the Xjp^ term 
in Eq. (5), which could be interpreted as the ‘pomeronic 
component’ of the fracture function. This contribution 
is clearly dominant for large /3, but fails to do so in the 
low-P and large-x^p region. 

In the analysis of both HI data setsi’i, which were 
obtained within different ranges of the variable t, we 
have assumed the universal validity of the exponenti^ 
behaviour in t measured by the ZEUS collaboration □ 
and have rescaled the data to a common range, given by 
that of the diffractive data of HI {—Mpx‘^/{1 — Xjp) = 
-tmin < -t<l GeV^). 

In Figure 3, we show the light-quark singlet and 
gluon fracture densities at the initial scale Qq and at 
an intermediate value of = 10 GeV^ for two char¬ 
acteristic values of x^p. The first one, Xjp = 0.005, 
which corresponds to well inside the diffractive regime. 


3 






































HI 1994 



0.2 0.4 0.6 0.8 0.2 0.4 0.6 0.8 

P 


Figure 3: Fracture function densities: a) and b) light-quark 
singlet, c) and d) gluons. 

shows a rather hard behaviour, whereas the second one, 
x,p = 0.2, which belongs to the leading-proton regime, is 
much softer. It is worth noticing that, at large the 
large-/? behaviour of the distributions is not well con¬ 
strained because leading-proton data are only available 
in the very small-/? range. 

As it is also shown in Figure 3, gluons carry much 
more impulse than quarks, specially in the case of the 
diffractive regime and at low scales. The evolution damps 
down this gluon dominance, suggesting a large fraction 
of valence-like gluons in the ‘pomeron’, which is not so 
apparent for ‘reggeons’. 

The main conclusion that can be drawn about the 
gluon density from the analysis is that the distribution is 
important at large /?, at variance with the one for inclu¬ 
sive structure functions, but its behaviour cannot be pre¬ 
cisely determined yet from the available data, in particu¬ 
lar, the exponent of the (1 —/?) factor in the parametriza- 
tion. 

4 Scale dependence 

As it has been said, the fracture function approach leads 
to very dehnite predictions about the scale dependence of 
the semi-inclusive cross sections. The rigorous factoriza¬ 
tion of the cross sections achieved within this formalism 
allows a precise perturbative QCD analysis of the scale 
dependence of the data, as is usually done for ordinary 
structure functions. 

In figures 4 and 5 we compare HI diffractive and 
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Figure 4: Scale dependence of HI diffractive data and the one 
obtained evolving the Set A parametrization. 

leading-proton data, at hxed values of Xjp and /?, as a 
function of Q^, with the evolved fracture function. The 
scale dependence induced by the evolution equations is 
perfectly consistent with the data. 

It is also worth noticing that within the measured 
range, the scale dependence is dominated by the homoge¬ 
neous evolution. The effects of the inhomogeneous term 

HI 1995 (Preliminary) 
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Figure 5: Scale dependence of HI leading proton data and the 
one obtained evolving the Set A parametrization. 
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become important only for low values of Xl where frag¬ 
mentation functions are larger. These_£jTects are, how¬ 
ever, beyond the range of present dataE3. 

5 ZEUS measurements 

In order to check the distributions obtained in the previ¬ 
ous section,in the following we compare results obtained 
using our best parametrization (Set A) with data sets 
presented by the ZEUS collaboration, which had not 
been included in the fit. 




Figure 6: The fraction of DIS events with a leading-proton for 
different ranges of t as measured by ZEUS 


ZEUS 1994 



Figure 7: ZEUS diffractive data, against the expectation 
coming from the fracture function parametrization. 


The ZEUS collaboration has measured the fraction 
of DIS events with a leading-proton in the final stataa. 
After the adequate rescaling of the parametrization for 
the fracture function at an average value of = 10 
GeV^ we obtain a remarkable agreement with the data 
in the common x^-range, as shown in Figures 6a and 6b. 
Notice also the fact that the parametrization interpolates 
fairly well the diffractive and leading-proton behaviour 
between HI kinematical regions (thick lines) where the 
data used in the fit belong. IruFigure 7 we compare 
ZEUS diffractive measurements □ with the outcome of 
the parametrization after the appropriate evolution to 
the mean scale value of the data = 8GeV^, and the 
already mentioned rescaling in t. 

6 Conclusions 

We have shown that an approach based on fracture 
functions motivates and allows a unified description of 
both leading-proton and diffractive deep inelastic cross 
sections. A simple parametric form for this function 
gives a very accurate description of the data available 
at present providing a smooth interpolation between the 
distinctive behaviours of the two regimes, also in accor¬ 
dance with ZEUS data. The analysis also hints at some 
non-perturbative features, such as a strong-gluon dom¬ 
inance in the ‘pomeronic’ component with a character¬ 
istic valence-like behaviour. Finally, our results verify 
that the scale dependence of the data agrees with the 
one predicted by the fracture function formalism. 
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